Study of the Cabibbo-suppressed decay modes 

D° 7T-7T+ and D° 

The FOCUS Collaboration* 

J. M. Link a , M. Reyes a , P. M. Yager a , J. C. Anjos b , 
I. Bediaga b , C. G6bel b , J. Magnin b , A. Massafferri b , 
J. M. de Miranda 13 , I. M. Pepe b , A. C. dos Reis b , S. Carrillo c , 
E. Casimiro c , E. Cuautle c , A. Sanchez-Hernandez , C. Uribe c , 
F. Vazquez , L. Agostino d , L. Cinquini d , J. P. Cumalat d , 

B. O'Reilly d , J. E. Ramirez d , I. Segoni d , M. Wahl d , 
J. N. Butler 6 , H. W. K. Cheung 6 , G. Chiodini 6 , I. Gaines 6 , 

P. H. Garbincius 6 , L. A. Garren 6 , E. Gottschalk e , 
P. H. Kasper 6 , A. E. Kreymer 6 , R. Kutschke 6 , L. Benussi f , 
S. Bianco f , F. L. Fabbri f , A. Zallo f , C. Cawmeld g , 
D. Y. Kim g , A. Rahimi g , J. Wiss g , R. Gardner 11 , 
A. Kryemadhi 11 , C. H. Chang 1 , Y. S. Chung 1 , J. S. Kang 1 , 
B. R. Ko 1 , J. W. Kwak 1 , K. B. Lee 1 , K. Chc-j, H. ParkJ, 
G. Alimonti k , S. Barberis k , M. Boschini k , A. Cerutti k , 
P. D'Angelo k , M. DiCorato k , P. Dini k , L. Edera k , S. Erba k , 
M. Giammarchi k , P. Inzani k , F. Leveraro k , S. Malvezzi k , 
D. Menasce k , M. Mezzadri k , L. Milazzo k , L. Moroni k , 
D. Pedrini k , C. Pontoglio k , F. Prelz k , M. Rovere k , S. Sala k , 
T. F. Davenport III*, V. Arena 111 , G. Boca m , G. Bonomi 111 , 

G. Gianini 111 , G. Liguori m , M. M. Merlo m , D. Pantea m , 
D. Lopes Pegna m , S. P. Ratti m , C. Riccardi m , P. Vitulo m , 

H. Hernandez 11 , A. M. Lopez 11 , E. Luiggi n , H. Mendez 11 , 
A. Paris 11 , J. Quinones 11 , W. Xiong n , Y. Zhang 11 , 

J. R. Wilson , T. Handler?, R. Mitchell 13 , D. Engh q , 
M. Hosack q , W. E. Johns q , M. Nehring q , P. D. Sheldon q , 
K. Stenson q , E. W. Vaandering q , M. Webster q , M. SheafF 

a University of California, Davis, CA 95616 
h Centro Brasileiro de Pesquisas Fisicas, Rio de Janeiro, RJ, Brasil 
C CINVESTAV, 07000 Mexico City, DF, Mexico 



Preprint submitted to Elsevier Science 



5 February 2008 



d University of Colorado, Boulder, CO 80309 
e Fermi National Accelerator Laboratory, Batavia, IL 60510 
i Laboratori Nazionali di Frascati dell'INFN, Frascati, Italy 1-00044 
g University of Illinois, Urbana- Champaign, IL 61801 
^Indiana University, Bloomington, IN 47405 
1 Korea University, Seoul, Korea 136-701 
3 Kyungpook National University, Taegu, Korea 702-701 
^INFN and University of Milano, Milano, Italy 
^University of North Carolina, Asheville, NC 28804 
m Dipartimento di Fisica Nucleare e Teorica and INFN, Pavia, Italy 
11 University of Puerto Rico, Mayaguez, PR 00681 
° University of South Carolina, Columbia, SC 29208 
p University of Tennessee, Knoxville, TN 37996 
q Vanderbilt University, Nashville, TN 37235 
1 University of Wisconsin, Madison, WI 53706 



Abstract 

Using data from the FOCUS (E831) experiment at Fermilab, we present a new 
measurement for the branching ratios of the Cabibbo-suppressed decay modes D° — ► 
7r~7r + and D° — > K~K + . We measured: 

T(D° -► K-K + )/F(D° -► 7T-7T+) = 2.81 ± O.lO(stat) ± 0.06(syst), 

T(D° -► K-K+)/T(D° -► K--K+) = 0.0993 ± 0.0014(stat) ± 0.0014(syst), and 

F(D° -► vr-vr+)/r(Z) -» K-vr+) = 0.0353 ± 0.0012(stat) ± 0.0006(syst). 

These values have been combined with other experimental data to extract the 
ratios of isospin amplitudes and the phase shifts for the D — ► KK and D — > tttt 
decay channels. 



1. Introduction 

In recent years, hadronic decays of charm mesons in many decay modes have 
been extensively studied. The theoretical models that have been developed 
mainly describe the 2-body decay modes and these models have led to sev- 
eral successful predictions. However the branching ratio ^poZ^-^^ nas > for 
a long time, been a puzzle of charm physics. The decay modes D° — > K~ K + 
and D° — > 7r~7r + are both Cabibbo-suppressed and in first order perturba- 
tive calculation both receive contributions from the same diagrams (external 

* See http://www-focus.fnal.gov/authors.html for additional author information. 
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spectator and exchange). To first order in the SU(3) flavour symmetry limit, 
the above branching ratio should be one. However this ratio is reduced by 
a factor 0.86 due to a phase space difference and increased by a factor of 
Uk/ fit) 2 = 1-49 because of the different decay constants of the kaon and 
the pion. An overall ratio of 1.29 is thus expected. Including SU(3) breaking 
effects, the expected ratio can increase to an upper limit of about 1.4 [1,2]. 

However, the measured ratio is close to 2.5 [3]. Both penguin diagrams and 
final state interactions (FSI) have been proposed as sources of such a high 
value. However the penguin interference [4] seems to be unable to explain this 
large value. Some theoretical models [5,6] propose the FSI as the solution of 
this puzzle. 

In this paper we present a new measurement of the T ^d^k-^+) branching 
ratio obtained using data from the FOCUS experiment as well as an isospin 
analysis of the D — > KK and D — > tttt decay channels. 

FOCUS is a charm photoproduction experiment [7] which collected data dur- 
ing the 1996-97 fixed target run at Fermilab. Electron and positron beams 
(with typically 300 GeV endpoint energy) obtained from the 800 GeV Teva- 
tron proton beam produce, by means of bremsstrahlung, a photon beam which 
interacts with a segmented BeO target. The mean photon energy for triggered 
events is ~ 180 GeV. A system of three multicell threshold Cerenkov coun- 
ters performs the charged particle identification, separating kaons from pions 
up to 60 GeV/c of momentum. Two systems of silicon microvertex detectors 
are used to track particles: the first system consists of 4 planes of microstrips 
interleaved with the experimental target [8] and the second system consists 
of 12 planes of microstrips located downstream of the target. These detec- 
tors provide high resolution in the transverse plane (approximately 9 /jm), 
allowing the identification and separation of charm primary (production) and 
secondary (decay) vertices. The charged particle momentum is determined by 
measuring their deflections in two magnets of opposite polarity through five 
stations of multiwire proportional chambers. 

2. Analysis of D° -> 7^+ and D° -> K~K + 

The final states are selected using a candidate driven vertex algorithm [7]. A 
secondary vertex is formed from the two candidate tracks. The momentum 
of the resultant D° candidate is used as a seed track to intersect the other 
reconstructed tracks and to search for a primary vertex. The confidence levels 
of both vertices are required to be greater than 1%. Two estimators of the 
relative isolation of these vertices are returned by the algorithm: the first 
estimator (Isol) being the confidence level that tracks forming the secondary 
vertex might come from the primary vertex, while the second estimator (Iso2) 



3 



is the confidence level that other tracks in the event might be associated with 
the secondary vertex. Once the production and decay vertices are determined, 
the distance L between them and its error <r L are computed. The quantity 
L / a L is an unbiased measure of the significance of detachment between the 
primary and secondary vertices. These variables provide a good measure of 
the topological configuration of the event, so that appropriate cuts on them 
reject the combinatorial background effectively. 

In addition to large combinatorial backgrounds, D° — > 7r~7r + and D° — > 
K~ K + decays are difficult to isolate because of a large reflection from the 
Cabibbo favored D° — > K~tt + decays. Particle identification requirements for 
each decay mode have been chosen to optimize signal quality. To minimize sys- 
tematic errors on the measurements of the branching ratios, we use identical 
vertex cuts on the signal and normalizing modes. The only difference in the 
selection criteria among different decay modes lies in the particle identification 
cuts. 

In the D° — > it~it + and D° — > K~ K + analysis, we require L / a L > 10, 
Isol < 10% and Iso2 < 0.5 %. We also require the D° momentum to be in 
the range 25 — > 250 GeV/c (a very loose cut) and the primary vertex to be 
formed with at least two reconstructed tracks in addition to the seed track. 
The Cerenkov identification cuts used in FOCUS are based on likelihood ratios 
between the various stable particle identification hypotheses. These likelihoods 
are computed for a given track from the observed firing response (on or off) 
of all the cells that are within the track's (/3 = 1) Cerenkov cone for each 
of our three Cerenkov counters. The product of all firing probabilities for all 
the cells within the three Cerenkov cones produces a x 2 Tike variable Wi = 
— 2 ln(Likelihood) where % ranges over the electron, pion, kaon and proton 
hypotheses [9]. All the kaon tracks are required to have = W n — Wk 
(kaonicity) greater than 1, whereas all the pion tracks are required to have 
= Wk — W v (pionicity) exceeding 1.5. Using the set of selection cuts just 
described, we obtain the invariant mass distributions for 7r~7r + , K~K + and 
K~n + shown in Fig. 1. 

In Fig. la the mass plot shows a broad peak to the left of the signal peak 

due to surviving contamination from D° — > K~ir + events. The shape of the 
reflection peak has been determined by generating Monte Carlo D° — > K~tt + 
events and reconstructing them as 7r~7r + . The mass plot is fit with a func- 
tion that includes a Gaussian for the signal, a third-order polynomial for the 
combinatorial background and a shape for the reflection background obtained 
from the Monte Carlo simulation. The amplitude of the reflection peak is a fit 
parameter while its shape was fixed. The low mass region is excluded in the fit 
to avoid possible contamination due to other charm hadronic decays involving 
an additional 7r°. A least squares fit gives a signal of 3453 ± lll7r~7r + events. 
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The K~K + mass plot, shown in Fig. lb, is fit with a function similar to that 
for the 7r~7r + fit. In the K~K + case, the reflection peak is on the right of the 
signal. The fit gives a signal of 10830 ± 148K~K + events. 

The large statistics K~n + mass plot of Fig. lc is fit with two Gaussians 1 with 
the same mean but different sigmas to take into account the different resolution 
in momentum of our spectrometer [7] plus a second-order polynomial. The fit 
gives a signal of 105030 ± 372K~n + events. 

The fitted D° masses are in good agreement with the world average [3] and 
the widths are in good agreement with those of our Monte Carlo simulation. 



3. Relative Branching Ratios 

The evaluation of relative branching ratios requires yields from the fits to 
be corrected for detection efficiencies, which differ among the various decay 
modes because of differences in both spectrometer acceptance (due to different 
Q values for the various decay modes) and Cerenkov identification efficiency 

From the Monte Carlo simulations, we compute the relative efficiencies to 
be: Jfefgy = 0.897 ± 0.003, ff^gg = 0.963 ± 0.003 and ffe^y = 
1.074 ± 0.004. Using the previous results, we obtain the following values for 
the branching ratios: ^ffp^ = 2.81 ± 0.10, PpCav) = °- 0993 ± 0-0014, 
and fId^kH) = °- 0353 ± 0-0012. 

Systematic uncertainties on branching ratio measurements can come from dif- 
ferent sources. We determine three independent contributions to the system- 
atic uncertainty: the split sample component, the fit variant component, and 
the limited statistics of the Monte Carlo. 

The split sample component takes into account the systematics introduced 
by a residual difference between data and Monte Carlo, due to a possible 
mismatch in the reproduction of the D° momentum and the changing exper- 
imental conditions of the spectrometer during data collection. This compo- 
nent has been determined by splitting data into four independent subsamples, 
according to the D° momentum range (high and low momentum) and the 
configuration of the vertex detector, that is, before and after the insertion 
of an upstream silicon system. A technique, employed in FOCUS and in the 
predecessor experiment E687, modeled after the S-factor method from the 
Particle Data Group [3], was used to try to separate true systematic varia- 
tions from statistical fluctuations. The branching ratio is evaluated for each 
of the 4 (= 2 2 ) statistically independent subsamples and a scaled variance a 



1 With the lower statistics of the K K + and tt tt + signals a single Gaussian gives 
an adequate fit. 
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(that is, where the errors are boosted when x 2 /(iV — 1) > 1) is calculated. The 
split sample variance Cgpi^ is defined as the difference between the reported 
statistical variance and the scaled variance, if the scaled variance exceeds the 
statistical variance: 



°split 
a split 



a 2 



stat 



if 
if 



° > CT stat 
d < a stat 



Another possible source of systematic uncertainty is the fit variant. This com- 
ponent is computed by varying, in a resonable manner, the fitting conditions 
on the whole data set. In our study, we changed the background parametriza- 
tion (varying the degree of the polynomial), the fit function for the reflection 
peak (the reflection shape from the Monte Carlo was replaced by a Gaus- 
sian), and the use of two Gaussian for the fit of the peak of D° — > 7r~7r + 
and D° — > K~K + . The values obtained by the various fits are all a priori 
likely, therefore this uncertainty can be estimated by the simple average of 
the measures of the fit variants: 



CT fit 



i=i 



xt 



N(x)< 



N — 1 



(2) 



Finally, there is a further contribution due to the limited statistics of the Monte 
Carlo simulation used to determine the efficiencies. Adding in quadrature the 
three components, we get the final systematic errors summarized in Table 1: 



Table 1 



Source 


a BR(KK/Kw) 


<7 BR(7nr/A'7r) 


Split sample 


0.0005 


0.0000 


Fit variant 


0.0012 


0.0006 


MC statistics 


0.0003 


0.0001 


Total systematic 


0.0014 


0.0006 



Sources of uncertainty on the p^-^+j* and branching ratios. 

The final results are shown in Table 2 along with a comparison with the 
previous determinations. 



r(7T-7T+) 



4. Isospin analysis of D — > KK and D — > nn channels 

Final State Interactions (FSI) can dramatically modify the observed decay 
rates and complicate the comparison of the experimental data with the the- 
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Experiment 


T{D°-*K-K+) 

r(D°^K-n+) 


r(D°^7r-7r+) 

r(D°^K-TT+) 


r(D°->K-K+) 

r(£> U -»7I— 7T+) 


E687 [10J 
E791 [11] 
CLEO [12] 
E831 (this result) 


0.109 ± 0.007 ± 0.009 
0.109 ±0.003 ±0.003 
0.1040 ±0.0033 ±0.0027 
0.0993 ±0.0014 ±0.0014 


0.043 ± 0.007 ± 0.003 
0.040 ± 0.002 ± 0.003 
0.0351 ± 0.0016 ± 0.0017 
0.0353 ± 0.0012 ± 0.0006 


2.53 ± 0.46 ± 0.19 
2.75 ±0.15 ±0.16 
2.96 ±0.16 ±0.15 
2.81 ±0.10 ±0.06 



Table 2 

Comparison with other experiments. 



oretical predictions. By means of the isospin analysis of the decay channels 
D — > KK and D — > tttt, we can gain some insight on the elastic component 
of the FSI (pure rotation in isospin space). 

Let us consider the D — > 7T7t transistions: -D — > 7r~7r + , D° — > 7r°7r° and D + — > 
7T + 7r°. The decay amplitudes can be expressed in terms of isospin If = (Ao) 
and If — 2 (A 2 ) amplitudes. The final state with isospin If — 1 is forbidden by 
Bose statistics for an angular momentum zero system of two pions. We denote 
by A + ~, A 00 and A +0 the decay amplitudes for the D° -> 7^+, D° -> 7r°7r° 
and -D + — > 7r + 7r , respectively. Expressing the decay amplitude in terms of 
isospin amplitudes, we have [13,14,15]: 



^ + ~ = ±\/^|A)|exp(^o) + \ -\A 2 \exp(i5 2 ) (3) 



i|A |exp (iS ) + J-\ 



A 00 = -\ l\Ao\ exp (iS ) + \ l\A 2 \ exp (i5 2 ) (4) 



A +0 = ±y||A 2 |exp^ 2 ) . 



(5) 



Adding the decay amplitudes in quadrature, we find the ratio of the magni- 
tude of isospin amplitudes and their relative phase shift difference in terms of 
measured branching fractions: 2 



M |2 
A 1 



2 
3 


| A 0+ 


2 


1 A+- 


2 + 1 


A 00 


2 


2 
3 


A+° 


2 



(6) 



, s r , 3 | A+- | 2 -6 | A 00 | 2 ±2 | A+° | 2 

cos(5 2 - S ) = — , . (7) 

4^2 | A+° P^fd A 00 | 2 + | A+- | 2 )- | A+° | 2 



2 The relationship between the isospin amplitude and the branching fraction is 
r u = g^r jp- | A l i | 2 , where is the center of mass 3- momentum of each final par- 
ticle [3]. 
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The decay rate T(D° — > 7r~7r + ) has been determined from our measurement 
of the branching ratio jv^~+) ; whereas the other D — > 7T7t decay rates and 
lifetimes have been taken from the Particle Data Group compilation [3] . 

The results are shown in Table 3. In contrast to K — > irn decays, where the 
transitions are dominated by the I = 1/2 amplitude (A I = 1/2 rule), the A 2 
amplitude in D — > 7T7t is comparable to the A amplitude. Furthermore, there 
is a large phase shift difference between the isospin amplitudes. According to 
Watson's theorem [16], this phase shift cannnot arise from the weak processes 
alone and thus constitutes direct evidence for FSI [17]. 

In the same way, we can consider the two-body D — > KK transitions: D° — > 
K + K~, D° -> K°K° and D+ -> if+K . The decay amplitudes A« for the 
£) — > K l Ki decay modes can be expressed in terms of A and A 1 isospin 
amplitudes [18]: 



A + =-^(|Ai|exp(i(Ji) + |i4o|exp(i<Jo)) (8) 
A 00 = i= exp {id,) - \A \ exp (i5 )) (9) 
A +0 = v^|Ai|exp(i(Ji) . (10) 



Using the previous decompositions, we can express the ratio of the magni- 
tudes of the isospin amplitudes and their phase shift difference in terms of the 
measured branching fractions: 



M |2 
A 1 







> 


2 | 


A+- 


| 2 +2 


| A 00 

A+- 


2 _ 


1 A+° | 
A 00 2 


2 




1 A+° | 


V 2 


| A 00 | 


2 +2 


| A+~ 


2 _ j ^+0 2 



cos(5i - 6q) = . / = • ( 12 ) 



The T(D° — > K~K + ) decay rate has been determined from our measurement 
of the branching ratio w^f+j ' ^ ne T(-D + — > A' + K°) from a previous mea- 
surement of FOCUS [19] and the remaining decay rate from the Particle Data 
Group compilation [3]. 

The results are shown in Table 3. Analogously to the D — > nn case, the two 
D — > KK isospin amplitudes are of the same order of magnitude, although 
the isospin phase shift difference is smaller. 

The isospin analysis of the D — > KK and D — > nn decay channels is sum- 
marized in Table 3 (the quoted errors are obtained adding in quadrature the 
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statistical and systematic errors) along with a comparison to previous deter- 
minations by CLEO [15,18]: 



Quantity 


CLEO 


E831(this result) 


1 M | / | A | 


0.72 ±0.13 ±0.11 


0.65 ±0.14 


<5 2 - 5 


(82.0 ±7.5 ±5.2)° 


(83.6 ± 10.0)° 


\A 1 \/\A \ 


o-6i±g:iS 


0.56 ± 0.04 


Si - S 


(28.4+^)° 


(37.1 ±7.5)° 



Table 3 

Isospin analysis for D — > KK and D — ► 7T7T decay modes, where | A2 | / | A) I and 
^2 — So refer to D — > tttt, while | | / | A | and Si — Sq to D — ► ifif . 

These results show that strong interactions, acting on the final particles, play a 
very important role in D — > KK and D — > 7T7T decays, modifying the measured 

r(7T 7T+) 

Another way to see the elastic FSI effect on the y^F^) branching ratio is to 
compute the ratio of the sums over the D° isospin rotated decay modes [20]: 
r [5 K t\ + V} K Q i Q ( ? ■ As opposed to the _ K 1) ratio, this ratio is not affected 

1 (IT 7T + ) + l (7T U 7T U ) rr 1 (IT 7T+) 

by elastic FSI. 

Using these measurements for ^(K-l+l ano - t(k-C+) ano - ^he PDG [3] values 
for the other modes, we compute: 



T(K-K+) + T(K°K°) 

r(7T-7T + ) + r(7T 7r°) 



= 2.06 ±0.24 



(13) 



This ratio is lower than the 



T(K~K+) 



branching ratio, but still above the ex- 



r(7r-7r+) 

pected value of 1.4. Therefore, the elastic FSI cannot account for all the dis- 
crepancy between theory and experiments. An inelastic FSI that also allows 
the transition KK — > ixix seems to be the most reasonable explanation [5]. 



4. Conclusions 

We have measured the following branching ratios: r(D°^ir-^+) ' r(D°^K-C+) 

an d T r(D0-^-n+) ■ comparison with previous determinations has been shown 
in Table 2. Our results improve significantly the accuracy of these measure- 
ments. 
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An isospin analysis of the decay channels D — > KK and D — > 7T7t shows that 
final state interactions play an important role in these hadronic decay modes. 
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Fig. 1. Invariant mass distribution for 7r _ 7r + (a), K~K + (b) and K~ir + (c). The fit 
(solid curve) for the Cabibbo-suppressed decay modes of D° is to a gaussian over a 
polynomial (for the combinatorial background) and a function obtained with Monte 
Carlo simulations for the reflection peak. 



12 



